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We investigate the localization of four-dimensional metastable gravity in supersymmetric dilatonic 
domain walls through massive modes by considering several scenarios in the model. We compute 
corrections to the Newtonian potential for small and long distances compared with a crossover scale 
given in terms of the dilatonic coupling. 4D gravity behavior is developed on the brane for distance 
very much below the crossover scale, while for distance much larger, the 5D gravity is recovered. 
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One of the main interesting results in higher dimensional physics in recent years, is the possibility of localization of 
\ gravity on braneworlds which was initiated in the seminal paper by Randall- Sundrum (RS) This attracted great 

■ attention because it was the first example where four-dimensional gravity could be localized even in the presence of 
£N| 1 extra dimensions with infinite size. Until then the well-known Kaluza-Klein mechanism to make compactification of 

extra dimensions with very small radii was the only largely considered alternative to achieve four-dimensional gravity 
from higher-dimensional theories. Braneworlds attached to extra dimensions with infinite size is able to localize four- 
dimensional gravity because the 3-brane is embedded in a curved higher dimensional space. More specifically, in the 
original set up [l[ , a 3-brane with tension sourced by a delta function is embedded in a warped five-dimensional bulk 
\ with negative cosmological constant known as AdS§ space. The relation between the brane tension and the bulk 
cosmological constant plays the role of producing a zero mode solution to the linearized Einstein equations. This 
solution given in terms of the metric developing the suitable fall-off is responsible for localizing the four-dimensional 
gravity. The following massive modes of the full spectrum also enter in the computation of the Newtonian potential 
describing corrections to the four-dimensional gravity. They should be highly suppressed at large distance but should 
deviate the Newton's law from the four-dimensional behavior at very small distance, implying that the five-dimensional 
C*~) ■ character of the full space shows up at very high energy. 

Since the RS paper, many extensions to thick 3-brane have appeared in the literature Q as domain wall solutions 
in five or higher dimensional theories with scalar fields coupled to gravity. Most of them are conceived in such a way 
that can be thought of as a bosonic sector of a supergravity theory or at least the sector of the sometimes referred 
to as fake supergravity 0] . One of the major problems in realizing the RS set up in five-dimensional supergravity is 
the difficulty to find normalizable zero modes. This issue was explored several times in the literature — see Q for a 
review. The main point is that the warp factor found does not enjoy the required fall-off along the extra dimension. 
!>• ■ One can rephrase this problem in terms of the volume of the five-dimensional spacetime. In the RS scenario, in spite 
KJ{ , of the fact the extra dimension is infinite, the volume of the AdS§ bulk is finite. In five-dimensional supergravity is 

■ hard to find solutions that gives finite volume. 

Surprisingly enough, soon after the appearance of the RS scenario, it was presented alternative scenarios by Gregory- 
Rubakov-Sibiryakov (GRS) [U an d Dvali-Gabadadze-Porrati (DGP) Q where the volume of the five-dimensional 
spacetime is not necessarily finite to localize gravity. This opened the possibility of localizing four-dimensional gravity 
in braneworlds embedded in asymptotically flat five-dimensional spaces. The price to pay is instead of zero modes 
the localization is due to massive graviton modes and the localization is metastable. In spite its metastability, the 
localization of metastable gravity can survive long enough if one probes distances very much small than a crossover 
scale. In these alternative scenarios, the brane is still sourced by a delta function but the bulk in both cases is 
asymptotically five-dimensional Minkowski space. 

Furthermore, one has also considered thick brane realizations of metastable localization of gravity in the literature 
" l5j . In most of these cases one has considered five-dimensional theories of scalar fields coupled to gravity in the 
realm of 'fake supergravity' where several issues were discussed. One should, though, also consider such realizations 
in supergravity domain wall solutions [l6j . 

In the present study, we shall consider the supersymmetric dilatonic solution that can be found from a higher 
dimensional supergravity theory that after specific compactifications , such as toroidal compactification, turns to a 
simpler five-dimensional theory of a scalar field coupled to gravity |17H2Cl | in the same fashion of those conceived in 
■ We shall investigate the possibility of localizing four-dimensional metastable gravity in supersymmetric dilatonic 
domain walls through massive modes by considering several scenarios in the model. 
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II. GRAVITY ON DILATONIC DOMAIN WALLS 



Let us briefly introduce the mechanism behind gravity localization in braneworlds. We consider a general dilatonic 
(D — 2)— brane solution in a spacetime in D dimensions discussed in Refs. [l7l-l20|. The spacetime is assumed to be 
the direct product Mn-i x K of the (D — 1)— dimensional spacetime Md_i and some noncompact internal space 
K. We shall focus our attention to obtain solutions describing braneworlds as the worldvolumc of domain walls. The 
Einstein-frame action for the dilatonic domain wall is 
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where kd is the D— dimensional gravitational constant. The cosmological constant term is multiplied by the dilaton 
factor which depends on an arbitrary dilaton coupling parameter a in any D-dimensional spacetime. The dilatonic 
domain wall solution given in terms of a conformally flat form is 



G MN dx M dx N = C(z)[-dt 2 + dx\ + ... + d 2 D _ 2 + dz 2 }, 
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In order to study the spectrum of the graviton in the dilatonic domain wall background, we consider the metric 
describing the small fluctuation h lllJ {x p ,z) of the (D — 1)— dimensional Minkowski spacetime embedded into the 
conformally flat D— dimensional spacetime 



g^dx^dx" = C(z)[(r]p V + h^dx^dx" + dz 2 ], 



(5) 



where Ih^l «l. 



By using the linearized Einstein equations in the transverse traceless gauge, i.e., d^h^ = 0, and the metric fluctu- 
ation in the form h^{x p ,z) = hffl(x p )C~( D ~ 2 ' > / 4: ipm{ z )i where rj^d^dvhffl = m 2 /i| 1 ™' ) , we obtain the Schroedinger- 
like equation 



3—5- 4>m(z) + U(z) l/j m (z) = m 2 1pm(z), 

a z A 



with the potential 



U(z) 
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where the prime denotes derivatives with respect to z. Now, we shall discuss the properties of the KK modes for 
different values of A taking Q > 0. We explore the Schroedingcr-like potential that arises by considering the following 
ranges and study the corrections of the Newtonian potential that are derived from wave functions that emerge 
according to each choice of A. Furthermore, there are no tachyonic modes because the structure of this potential 
allows to write the Schroedinger-like equation in a quadratic form form [H 



dz 



Hip(z) = 

where A(z) — \ In C(z). The zero mode obeying H ipo = is given by 
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It is easy to conclude from the equation Q that for Q > there exists normalizable zero mode only in the range 
A < — 2, fl9l. [20|. In the following we shall focus on four cases that do not fall into this range such that there is no 
zero mode. This is to investigate how easy four-dimensional massive gravity can be localized on the brancworld. 
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III. FOUR-DIMENSIONAL GRAVITY ON DILATONIC DOMAIN WALLS 

In the present section we shall focus on four scenarios to four-dimensional gravity on dilatonic domain walls by 
considering four ranges of A. 

• i) A > case. 
In this case, let us use the potential in Eq. ([7]) given by 

nt \ Q 2 (A + 1) ,2Q 

U(z) = 9+— z-o(z), (10) 

A2(1 + A±2Q| Z |) 2 A 
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FIG. 1: The potential Q as a function of z, for A = Q = 1. 
whose general solution of the Schroedinger-like problem for massive modes is given by [2l| 



J n \m\ - + \z 



F(m)N a [m[- + \z\ 



(11) 



where J n and N n are Bessel functions of the first and second kind, respectively, n 



k 



- See 



4 (A+2)^ 1 

Fig. [T] In order to obtain the correction term to the four-dimensional Newton's law between two unit masses on the 
3-brane (D = 5), one needs to obtain the probability density of massive modes on the brane. The asymptotic behavior 
of \4>m (0) | 2 can be found from the argument in the Bessel functions, the normalization factor C m and F{m) given by 



\C m \ 2 = 



2 (1 + F 2 (m))' 



(12) 
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Now, we use the Lommel's formula, J n +i(x) N n (x) — N n+ i(x) J n (x) = 2/irx, to find the asymptotic forms of Bessel 
function as J n (x) ~ y^^ cos ( x ~ l/2n7r + l/4 7r) and iV n (x) ~ ^/^sin (a; — 1/2 n7r + 1/4 7r) for x 3> 1, J n (a;) ~ 
(l/2x) n and -/V n (a;) ~ — 1/n (l/2x) n for .t <C 1, such that the probability density assumes the simplified forms 

'^ (0)|2 ~ ^fr^) ' m>>fc ' (14) 

2 /m\2n+l 

1^ (0)| 2 ~ — -2 (-) , m«fc. (15) 
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To compute the correction to the four-dimensional Newtonian potential generated by the massive modes we use the 
now well-known formula [l[ 

V (r) = -l-\ik(0)\* + SV(r), (16) 
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where the first term is the contribution of the zero mode and the second term corresponds to the corrections generated 
by the exchange of massive Kaluza-Klein modes 

SV (r) = ^ dm\ip m (0)\ 2 e~ mr - (17) 
r Jo 

In our present study we have no zero modes, then we focus our analysis in two distinct regions from the integral in 
the second term only, where we integrate out the probability densities for the continuous massive modes on the brane 
at z = 0. We divide this integral into two regions limited by the crossover scale 

V(r) = SV(r) = ^ ( — ^ I" f^) 2n+1 e-mr dm+ [°° 2m * -mr dm \ (18) 

r \(e + l) 2 J \kJ 7i Trfc 2 ( e 2 + ^.) J> v ' 

where wc define the crossover scale r c = 1/k. Thus, the Newton's law correction for massive modes at a distance r 
is given as follows. In the limit of the crossover scale is very small, i.e., r c — >• (A <C Q), the first integral in (|18p is 
dominant and gives the following answer 



V(r) = 



(e + iy M 5 3 r* 



M- 1 



(19) 



which is defined in terms of Kummer's confluent hypergeometric functions M=M^n, n + ^ ; . We found the familiar 

behavior of Randall- Sundrum scenario pj for the critical values of the dilaton coupling a 2 ~ 4 rgrE^A m this limit, so 
that 

8A(A + 2)' /Q 

Q (5 A 2 + 20 A + 24) \M 5 3 r 3 J 

On the other hand, for the crossover scale being very large, i.e., r c — > oo (A 3> Q) the second integral in (|18p 
is dominant and we can get approximately the form of AV(r) at a distance r, using the relation Ei(ix) = ci(x) + 
i [1/2 7r + si(x)], where ci(x) = — J°° cos tjt and si(x) = — sin i/t, as follows 



( sin fe) Ci fe)~ C ° S fe) si (^ 



Vij) ~ r 2 M 5 3 7r + 2rr c nM§ ' (21) 

We are now in a position to make comparison of the distance r with the crossover scale r c in order to analyze the 
large and small distance behavior. We use the asymptotic forms: ci(x) ~ -y + ln(x) and si(x) ~ x — 7r/2 for x <C 1, 
ci(x) ~ sin(x)/x and si(x) ~ — cos(x)/x for x ^> 1. 

For small distance, i.e., r/r c -C 1 we obtain the following form 



F(r) -27^l-^ + 1 "l-JJ-- + 2 + °< ri )J- < 22 » 

where 7 « 0.577 is the Euler-Masceroni constant. Interestingly, at short distance the computed potential has the 
correct 4D Newton's law with 1 jr scaling. The next to leading correction is given by the logarithmic repulsion term 
in (p2|. 

On the other hand, for large distance, i.e., r/r c ^> 1, the potential in Eq. (|2"Tj) gives 

y W ~ 9 2 lr 3 ■ r » r - ( 23 ) 
2 r z M5 7r 

which describes the laws of 5.D gravity [fl 0] . 
• ii) — 1 < A < case. 

The following analysis is pretty similar to the previous one. In the present case, the Schroedinger-like potential ([7]) 
gives 

ru \ Q 2 ( A + !) 2Q„ . 
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FIG. 2: The potential as a function of z, for A = -1/2, and Q = 1. 



and the Schrocdingcr-like equation now gives the following solution for massive modes [21 



m (1 + a \z\) 



(25) 



in terms of the Bessel functions of first and second kind, Jp and Np respectively, where a = — - — 
— See Fig. [2] The jump condition at z = 
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(0) , 



(26) 



is obtained from the Schrocdingcr-like equation by using the properties of the delta function. The asymptotic behavior 
of \ip m (0) | 2 depends on the magnitude of the argument in the Bessel functions, the normalization factor and 
F 3 (m) given by 



IK 



( m ) |2 



2 a (1 + F 2 (m))' 



(27) 



, , ( -2 Js+i (s) m + e 3 Jb(-)) , , 

F 3 m = 4 \ a \\ , 28 

where e 3 = 2Q. Then, as in the previous case, using the Lommel's formula, J n+ i(x) N a (x) — N n+ i(x) J n (x) = 2/ttx 
and the asymptotic form of the Bessel functions, the probability density assumes the simplified forms 

, , ,„.,o 7T m 2 1 

W>™ ~ -t— ^ , 2 ,A 2V m»-; 29 

/ o „,2(l-,3) \ i 

fc(0)| 2 H U 2/?+ \ m<-. (30) 



3 (2a-e 3 ) 



Again, in this problem we have no zero mode, consequently, all the contribution to the Newtonian potential comes 
from the massive modes that we shall integrate out in the two regions below 

V(r) = SV(r) = ^- | _2 / m 2£+i e -mr dm + / m e- mr dm , (31) 

r \3 (2a-e 3 ) 2 Jo 4a 2 J_l (e 2 +4m 2 ) ^ 

where we define the crossover scale r c = 2/e 3 . In the limit of the crossover scale is very small, i.e., r c — > (A — > 
necessarily), the first integral in pip is dominant and we naturally find the familiar behavior of the Randall-Sundrum 
scenario [l[ 

V(r)~m^?. (32) 
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On the other hand, for the crossover scale being very large, i.e., r c — > oo (Q ^> 1) the second integral in (l3lT) is dominant 
and we can find AV(r) at a distance r, using the relation Ei(i x) = ci(x) +i [1/2 7r + si(x)], where ci(a;) = — J°° cos t/t 
and si(x) = — sin t/t, as follows 



4 Q 2 (2 + A) 2 I 4r 2 M 5 3 7r 



| e 3 (sin (y ) ci ) - cos (y) si (y)) 



16r7rA/| 



(33) 



We can now discuss the large and small distance behavior comparing with the crossover scale. 
For small distance, i.e., r/r c <C 1 (Q <C 1), we obtain the following form 



V(r) 



£3 



64 (2 + A) Ml r Q 2 



0(r 



(34) 



Notice that at short distances the potential has the correct AD Newtonian 1/r scaling, with a subsequently correction 
by the logarithmic repulsion term in (|34l) . 

Finally, for large distance, i.e., r/r c 3> 1, the potential in Eq. (|33[) gives 



v{r) + A" 



1 



1 

^2' r > 7 'c, 



(35) 



64 (2 + A) 2 Q 2 Vr 2 M 5 3 , 
which is in accordance with the laws of 5D gravity Q . 
• hi) A = — 1 case. 

The potential ([7]) now reduces to U(z) = —2QS(z) in this case — See Fig. [3] The probability density for the 



-1 >P- 

-1 



FIG. 3: The potential (0 as a function of z, for A = —1, and Q = 1. 

massive modes is given in terms of the scattering states governed by \ip m (0) | 2 that depends on the magnitude of the 
transmission T or reflection R coefficients. As usual, the jump condition at z = 0, 



dz 



tpm (z) 



= -2QiP (m) (0) 



(36) 



is obtained from the Schroedinger-like equation by using the properties of the delta function. We now consider the 
general wave functions for scattered states in the form 



c lKZ +Re 
Te lKZ , 



~ lK \ z<0 
z>0 



(37) 



where k is the wave number. The simplified form for the probability density in this case is 



\T\ 



hMo) r 



2A 



(38) 



Now we proceed as in the previous cases. Again, all the contribution to the Newton's law arises from the continuous 
massive modes on the branc such that we have 



V(r) = SV(r) 



Mr 



3 /■oo 
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r dm 1 



(39) 
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where we define the crossover scale r c = l/y/2A. Notice this asks a bulk positive cosmological constant. 
For the crossover scale being very large, i.e., r c — > oo, the integral in (|39j) turns 



V(r) 
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For small distance, i.e., r/r c <C 1 we obtain the following form 

1 



V(r) 



rr c M 5 A \r c 



r I I r 

7 + m 



+ \ + °( r ^ ) > r ^ Tc 



(40) 



(41) 



As we have found, in the previous cases, at short distances the potential has the correct AD Newtonian 1/r scaling. 
On the other hand, for large distance, i.e., r/r c S> 1, the potential in Eq. (j40]) gives 



V(r) 



1 



r > r c 



r 2 M 5 3 7T ' 

which naturally signalizes a 5D gravity behavior [H, Q. 

• iv) —2 < A < — 1 case. 
In this last case, the Schrocdingcr-like potential (JT]) is formaly the same as given in (pM)) . i.e. 



(42) 



U(z) 
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IA + 11 
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A 2 (l-(^)Q|,|) IAI 
which massive modes solution of the Schroedinger-like problem is given by [2l| 



(43) 



(z) = yiT^N (j* ( m(1+ ^ N) ) + *k(m) JV* 



m (1 + «i |z| 
ai 



(44) 



in terms of the Bcssel functions of first and second kind, Jp 1 and Np 1 respectively, where a.\ — ^ A | A ^ — , (3\ 

1 ^/8-8|A|+A2 



2 (|A|-2) — ' This problem is similar to that of a particle trapped into an infinite box up to a delta function that 
we use to impose boundary conditions — See Fig. [4] The jump condition at z = 0, is given by 
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FIG. 4: The potential |7| as a function of z, for A = -3/2, and Q = 1. 



z=0 



-m- Wm (0) 



(45) 



The asymptotic behavior of \ip m (0) | 2 are given in terms of the normalization factor Kg" 1 '' and Fi(m), i.e. 



IK 
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-2J fa+1 (^) m + 6 4 J gl (t)) 



(46) 
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where £4 = |2jr (— |A| + 2 /3i (|A| — 2)). Thus, the asymptotic probability density assumes the simplified forms 



'^ (0)|2 ^ 4^f(^Tw)' m>> ^ (47) 



> 2 1 



l^(0)| 2 ~ II ° ia(1 " ft L )™ 2 * +1 , m«I 



and the boundary conditions of the potential applied to the wavefunction, Eq. (j44|) , 

^ m )(—)=0, z>0; V (m) f -— J=0, z<0 (49) 



allows us to determine how the parameter £4 is related to the graviton masses, i.e. 

2 



2 1 2 / J/31 1^ 



4^ ^ 



V-W (fr) 



For qj -c 1 (or in terms of the bulk cosmological constant |A| 3> 1), in the asymptotic regime the solution of the 
Eq. (|50[) is approximately given by 



£4771 \ai 



m (2/3i + 1)tt 



—L, k = 1,2,3... . (51) 



In such a regime, it is reasonable to approximate the potential generated by discrete massive graviton states as a 
summation of Yukawadike potentials, which makes the total effective potential to have the form 0, 0, Q 

p -m r p-m r 

V(r) = ^(0)^— + £ tf£(0)-^, (52) 

n 

where the first term is contribution of zero mode and the second term corresponds to the correction term which 
is generated by the exchange of KK-modes. The set of discrete states obtained, however, may be replaced by a 
continuous treatment for ol\ <C 1, such that — > J. In this problem we have no zero modes, then our analysis focus 
on in two distinct regions from the second term. As in previous cases, we integrate out the massive modes as follows 

V { r)=6V(r) = OS™, rm^e— dm+ « f ^ e — dm), (53) 

r \3 (2 ai -£ 4 ) 2 Jo Aa iJh ( £ 4 + 4m ) / 

where we define the crossover scale r c = 2/e4. 

Once again, in the limit of the crossover scale is very small, i.e., r c — > (A — > —4 + 2\[2, — ^ 0), the first 

integral in (|53|) is dominant and we find the familiar behavior of Randall- Sundrum scenario [l[ 

6 r 3 Ms 

On the other hand, for the crossover scale being very large, i.e., r c — > 00 (|A| <C Q) the second integral in (|53l) is 
dominant and we obtain the potential as follows 

VM 1 " A2 f , ^(Bin(^)ci(^)-cos(^)si(^)) \ 

l) ~4g2( 2 + A) 2 ^4r 2 Af 5 3 / 16r7rM 5 3 J - 1 j 

For small distance, i.e., r/r c <C 1 we obtain the following form 

^>~ «4,, + W £(^ +h re))-^ + i + °'^- r<<r < (56) 

Notice that, at this limit the potential has the correct AD Newton's law with 1/r scaling. 
Finally, for large distance, i.e., r/r c ^> 1, the potential in Eq. ([5"5j) gives 

that recovers the laws of 5D gravity [!, @ . 
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IV. CONCLUSIONS 



In this paper, by considering massive graviton modes coming from a dilatonic domain wall solution of five- 
dimensional supergravity, we have shown that the gravitational potential corresponds to the usual Newton potential 
which scales with (— > 1/r) at short distance and has a five-dimensional behavior scaling with (— > 1/r 2 ) at large 
distance compared with the crossover scale which depends on the parameter A chosen in several ranges worked in the 
present analysis. This study showed that from a five-dimensional supergravity theory with a scalar field describing 
the dilaton, the emergence of four-dimensional gravity on a 3-branc is possible even if the brane is embedded in an 
asymptotic five-dimensional flat space, below a crossover scale and the manifestation of extra dimensions does not 
necessarily occur only at short distances as commonly expected. Localization of gravity in other dilatonic domain 
walls configurations has also been addressed in [pj [23| ■ 
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